A Note on Perturbations During a Regular Bounce 
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We point out an inconsistency in a method used in the literature for studying adiabatic scalar 
perturbations in a regular bouncing universe (in four dimensions). The method under scrutiny 
consists of splitting the Bardeen potential into two pieces with independent evolutions, in order to 
avoid a singular behavior at the boundaries of the region where the null energy condition (NEC) is 
violated. However, we argue that this method violates energy-momentum conservation. 

We then introduce a novel method which provides two independent solutions for the Bardeen 
potential around the boundaries, even in the case of adiabatic perturbations. The two solutions are 
well behaved and not divergent. 
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In recent years, regular toy-models of a bouncing uni- 
verse have received a lot of attention due to string in- 
spired singular bouncing models such as the cyclic sce- 
nario. Of particular interest is the evolution of scalar per- 
turbations, since their spectrum is directly observable in 
the anisotropies of the cosmic background radiation and 
the large scale structure of the universe. The majority of 
toy models are four dimensional, requiring at least two 
matter fields, of which at least one has to violate the null 
energy condition (NEC). 

If one focuses on adiabatic perturbations, the evolu- 
tion equation for the Bardeen potential becomes singular 
at the boundaries of the NEC violating region before and 
after the bounce. One way out of this dilemma was advo- 
cated in pj and subsequently used e.g. in 0- [J [!] : One 
splits the Bardeen potential in two components, each of 
which satisfies a regular second order differential equa- 
tion. 

In this note we first show that this method is incon- 
sistent because the fluid conservation equations are vi- 
olated. Therefore all models using this method have to 
be re-evaluated. We then introduce a novel method, pro- 
viding two independent and well behaved solutions for 
the Bardeen potential around the boundaries of the re- 
gion with NEC violation, even in the case of adiabatic 
perturbations. 

For simplicity, we work with a two fluid model with 
T v = T (a) v ± T (b) v and equations of state p {l) = wip (l) 
for I = a, b. Note that a bounce for the scale factor a 
occurs only in case of a negative sign in front of 
Perturbing the metric 



{b)v 



ds 2 = a 2 [(l + 2$)d?7 2 - (1 -2<$>)S ij dx i dx j ] , (1) 

where $ is the Bardeenpotential (longitudinal gauge, no 
anisotropic stress, see [g for details), and perturbing also 
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the energy momentum tensors 

far" ) - om ( ( 1 + w i) v d)^ \ (2) 

where Sp^ = Wi8p {l) and := 5p( l} /p {l} , the perturbed 
Einstein equations for /C = (a spatially flat universe) 
read 



V 2 <& - mm + §') = -K 2 ( P{a) t {a ) ± P(6)£(6)) , 



(3) 



±Wb£(b)P(b)) > ( 4 ) 



[*W + - --K 2 ( P{a) V {a)>i (l + W a ) 

±P(b)V { b),i( 1 + W b)) > 



(5) 



with k 2 = 8ir/M 2 . The energy conservation equations 
for each fluid, assuming no non-gravitational interactions 
between the fluids, are 

P(l) ((1 + Wl ) [V 2 V (l) - W] + = . (6) 

In the long wavelength limit the above relation yields 

3«> 



1 + wi 



3<J>', 



(7) 



which is another way of stating entropy conservation, 
that is 



S' 



(b) 



l+W a 1 + U'b 



0. 



(8) 



One can simply impose the adiabaticity condition by set- 
ting initially S(r]i n ) := and thus @ 



f(a) _ £(6) 



1 + w a 1 + Wb 



(9) 
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has to hold. 

Substituting this result back into and |0J and then 
combining the two, one arrives at a second order equation 
for <E>: 



= $" 



+3H$' 



Pa(l + W a ) ± + W b ) 



Pa(l + Wa) 2 ± Pb(l + Wbf 



(10) 



- {w a (w a + l)p a ± Wb(Wb + l)/9b) V 2 

+2H'( /9a (l + w a )± / 9 6 (l + u; b )) 

+n 2 ( Pa (l + w a ){l + 3w a ) ± p 6 (l + it; 6 )(l + 3w b )) 



4>. 



In case of a minus sign, which is needed for a bounce to 
occur, or some negative Wi this equation becomes singu- 
lar at the boundaries of the NEC violating region (see 
also Q). 

In the appendix of Q a split of <!> into <& a + <E>b is 
suggested, such that each satisfies 



(ii) 



+ 3H$; + (2H' + H 2 )<i> t = -k 2 w£ (1)P(1) .(12) 

By adding these equations one arrives at the original Ein- 
stein equations 10 and l@J). Combining l(TT|) and (|T2"|) one 

can derive for each the equation 

= &{ + 3H(l + wi)& l 

+ (-wiV 2 + 2H' + (1 + 3wi)H 2 ) . (13) 

Note that these equations are regular for each . This is 
already the first hint that the method is doubtful, since 
the singular behavior vanished miraculously. Another 
reason to doubt this method is the fact that by just us- 
ing equations © and Q), with no extra constraint like 
conservation of energy for each fluid or adiabaticity, one 
seems to be able to calculate the evolution of each 
and subsequently <I> itself. However, it is obvious that in 
the equations and Q three unknowns appear, so that 
two equations are insufficient to calculate their evolution. 
Therefore one can not expect the solutions of l|ll|) and 
(|12fl to be consistent with the conservation of the energy- 
momentum tensor or other constraint equations. 

One easy way to see this inconsistency is to look at 
the solutions in the long wavelength limit where equation 
(lllll simplifies to 



(14) 



Adding the time derivative of the above equation to TL x 
itself and 3Hx yields 

= 3#' (p (a )(l + u; )+p (6) (l + u>6)) , (15) 
P(6)£(6) = W b (p w (l + w b ) + P(a) (l+w a )) , (16) 



where we also used the background Einstein equations 
and p\ — — 3H(1 + WijPi. This result, together with the 
conservation equation JJJ, implies that the solutions of 
(|13fl have to satisfy 



$1 



P( a )(l+U> a ) P(b)(l+Wfc) 



(17) 



One might think that this constraint is satisfied in the 
case of adiabatic perturbations, but it is not as we shall 
see now. Using © in 1(14^1 yields 



n<s> a + K 



P(a) (1 + W a ) p(b) (1 + Wb) ' 

which can further be simplified by using l|17|) to 



P(a)(l + W a ) p( b )(l + Wb) 

This together with (|17fl implies 



a 



(18) 



(19) 



(20) 



yielding <& a tx However, this is clearly in contradic- 
tion to (|19l) . since the densities have a different depen- 
dency on conformal time if w a ^ Wb- Therefore the split- 
ting method itself is inconsistent and can not be trusted 
to regularize (jTU|l . 

We propose another mathematical technique that can 
be used instead to approximate the solutions of IjlOl) in 
the vicinity of p to t + Ptot — at ?7„ ec . Equation (|10|) is 
a second order differential equation that in fourier space 
(suppressing the subscript on <&k), has the following gen- 
eral form: 



A(r])<f>" + B(rj)<P' + C(k, 77)$ = 0, 
where A and B are related via 



(21) 



B = -A 1 , (22) 
since energy conservation for each fluid requires |10| 

f i l = -2Mp l {l + vn). (23) 

Our first goal is to derive one of the solutions around 
Vnec perturbatively. We can then obtain the second so- 
lution by means of the Wronskian method. 

Since a(j) nec ) ^ we can Taylor expand pi and conse- 
quently A, B and C around rj nec : 



A(6) = A^ + A 2 S 2 ^ s± 3 c 
B(S) = -Ai - 2A 2 S - 3A 3 6 2 
C(S) = C + C 1 5 + C 2 S 2 + .. 



A,S 3 



(24) 
(25) 
(26) 



where we defined S = r\ — rj nec and used A(r] nec ) = 0. 
We assume that A\ , corresponding to the linear term of 
A(5), does not vanish. This is usually the case and plays 
a crucial role in our analysis. 



3 



Furthermore, we assume the existence of an analytic 
solution around r\„ ec (we will see bellow that it is a jus- 
tified assumption |ll|). so that it can be written as 



$1 = ao + ot\8 + a 2 5 + 



(27) 



Substituting relations <E3I>- <E3> into (2U we obtain the 
following relations from the zeroth and first order equa- 
tions in 5 

0{6 ) = -Aia x +C a , (28) 

0{5 l ) => 0= {-2A 2 + C ) ai + da . (29) 

In general, this can only be satisfied if ao = ct\ = 0. 
Fortunately, this does not imply $i = 0, because the 
equations of higher order in 6 can all be satisfied recur- 
sively. In fact, we can compute a complete power series 
solution for $i 



0(S 2 ) 



a 3 = - 



2A 2 - Cp 
3A 1 



« 2 



(30) 



0(6 n ) => a n+1 -. 



S"=2[*( 3 + n~ 2i)A n+2 -i - C n ^i]ai 



(n- l)(n + l)Ai 



with the consequence 



2A 2 - Cp 
3^i ' 



(31) 



Knowing one of the solutions of 1|21|) . <&i, we can easily 
obtain the other solution, $ 2 , by using the Wronskian 
technique. The Wronskian for a second order differential 
equations is defined as 



(32) 



where $1 and $2 are the independent solutions of (|21ll . 
Henceforth, one can calculate $2 in terms of W and $1: 



W 

<!>j - <1>, / -^drj. 



(33) 



W itself satisfies the first order differential equation 

A{rj)W' + B{rf)W = . (34) 

By substituting B(rf) from ill'L'l) we can solve the above 
equation for W to 



W(t)) = PMv), 



(35) 



where (3 is just a constant. Combining this result for 
W(rf) with l|33(l and our solution for <I>i from H31(l . we 
end up with 

q 2 = _^_8A 2l -C 1s + A3§2 2 
2 

Note that although this solution is not analytic at 7/ = 
Vnec (5 = 0), it is well behaved in the sense that both, the 
solution and its first derivative, remain continuous and 
finite. Thus, we have constructed the approximate form 
of the two independent solutions of (|2*T)l or subsequently 
(|10|l . These can be used to match the solution on different 
sides of r] nec . 

To summarize, we have shown explicitly that the split- 
ting method, first introduced in [l| and used to show 
the regularity of the Bardeen potential, is intrinsically 
inconsistent. All models using this method have to be 
re-evaluated |T3 |. e.g. by using the technique introduced 
in this draft or by including entropy perturbations, as 
emphasized in Q and later on in (they derived a reg- 
ular forth order equation for the full Bardeen potential). 
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